Saari's homographic conjecture for planar 
equal-mass three-body problem in Newton gravity 



Toshiaki Fujiwara^, Hiroshi Fukuda^, Hiroshi Ozaki^ 
and Tetsuya Taniguchi'' 

^ ^College of Liberal Arts and Sciences, Kitasato University, 1-15-1 Kitasato, 
Sagamihara, Kanagawa 252-0329, Japan 

•^General Education Program Center, Tokai University, Shimizu Campus, 3-20-1, 
Orido, Shimizu, Shizuoka 424-8610, Japan 

E-mail: ^fujiwara@kitasato-u. ac.jp, ^fukuda@kitasato-u. ac.jp, 
ozaki@tokai-u.jp, *tetsuya@kitasato-u. ac.jp 

Abstract. Saari's homographic conjecture in A^-body problem under the 

Newton gravity is the foUowring; configurational measure fi = \/7U, which is 

the product of square root of the moment of inertia / = (XI '"fc)"^ 

and the potential function U = ^ miirij jrij, is constant if and only if the motion 

is homographic. Where represents mass of body k and r^j represents distance 

between bodies i and j. We prove this conjecture for planar equal-mass three-body 

problem. 

In this work, we use three sets of shape variables. In the first step, we 
use ^ = 3(j3/(2(g2 — 91 )) where q^, £ C represents position of body k. Using 
T\ = r2s/ri2 and r2 = r'3i/ri2 in intermediate step, we finally use /i itself and 
p = /'^/^/(ri2r23r3i). The shape variables and p make our proof simple. 
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1. Saari's homographic conjecture 

In 2005, Donald Saari formulated his conjecture in the following form [10l[Tl]; in the 
A''-body problem under the potential function 
\ ^ mim.j ^ 

l<i<j<N *J 

a motion has a constant configurational measure 

H = u (2) 

if and only if the motion is homographic. Here, r^j represents the mutual distance 
between the bodies i and j, and / represents the moment of inertia 

I = i ^f^y^ E ^i^j^fj. (3) 

l<*:<Af i<i<j<N 

Florin Diacu, Toshiaki Fujiwara, Ernesto Perez-Chavela and Manuele Santoprete 
called this conjecture the "Saari's homographic conjecture" and partly proved this 
conjecture for some cases [2] . Recently, the present authors proved this conjecture for 
planar equal-mass three-body problem for a = 2 [7. In this paper, we extends our 
proof to a = 1 , the Newton gravity. 

In section [21 we derive the equations of motion for the size change, rotation and 
shape change. To do this, we use the shape variable (, 

C = |^, (4) 

introduced by Richard Moeckel and Richard Montgomery [8]. Here, g/c G C, fc = 1, 2, 3 
represents position of the body k. Then, in the sectionlH we investigate motions with 
IJ, = constant and non-homographic, and we derive a necessary condition that must be 
satisfied by such motion. The contents in the sections[2]and[3]are review of our previous 
paper [3], although we changed few notations. To prove the Saari's conjecture, we will 
show that no finite orbit satisfies the necessary condition. To attain this purpose, the 
expression of the necessary condition by ^ is too complex. To simplify the expression, 
we will use other set of shape variables, 

ri = IC-1/21 =r23/ri2, = |C + 1/2| - r3i/ri2. (5) 

Then, using the invariance of the system under the permutations of {91, <Z2, ^a}, we 
rewrite the necessary condition in another set of shape variables fj, itself and p, 

\ri2 r23 T-ix) TxiTi'ir-ii 

that are manifestly invariant under the permutations. Since, we are considering [i, = 
constant orbits, variables fj, and p make our proof easy. This expression is given in 
section [H The proof of the Saari's conjecture is given in the section [5] In section [6j 
we give discussions. 
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2. Equations of motion 



In this section, we summarize ttie equations of motion for a = 1 in terms of size, 
rotation and shape. We don't assume fi = constant in this section. 

Let Qk & C he the position and mass = 1 for k = 1,2,3. We take the center 
of mass frame, '^qk =0. The Lagrangian is given by, 

2 



L 



2 ^ 



dqk 



dt 



U. 



We take the shape variable ( E C in l^. 
scaling and rotation, qk — ?> /^e^^qk with A, 6* G K 
us define — qk/{q2 ~ ^i)- Then, we have. 



(7) 

This variable is invariant under the 
Thus, C, depends only on shape. Let 



2 3' 



6 



1 C 



2C 



2 3'^^ 3 



(8) 



Since, the triangle 91 92 93 and ^i^2^3 are similar and have the same orientation, we 
have two variables / > and 6* S M, such that 

a- 



qk 



Vie 



We take /, 9 and C for dynamical variables. In the following, we identify C 
and X = {x,y) G M^. By direct calculations, we obtain the Lagrangian 



r I 

L = h - 

8/ 2 



|x A X 



Here, ' represents time derivative, x A x 
/x(x) = 



xy- 



1 2, ,, 

2 3'' 



1 



2(l + ||xr 
yx and 

1 



(9) 

+ iy 

(10) 



V/(X- 1/2)2 + y2 ^{X + 1/2Y 

Since, 6 is cyclic, the angular momentum C is constant of motion. 



.(11) 



de 



lie 



1 



ix A X 



Therefore, the total energy E is given by 



/2 C2 / 

= \ \ 

8/ 2/ 2 1 



(12) 



(13) 



The three terms in the kinetic energy are kinetic energy for the size change, for the 
rotation and for the shape change respectively. The equation of motion for / yields 
Lagrange- Jacobi identity, / = AE + 2U . From this equation, we get the following 



"Saari's relation" |10,, 



d_ 
di 



2 (1 



||X|2)2 



Using the 'time' variable s defined by 



ds 



1 



dt 2/ V ^ 3 
the Saari's relation is written as 



(14) 
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The equation of motion for x in terms of s is 

dx 



^^'V'^Ts (dy dx\ , ^^r-^^, 



(17) 



ds^ l + ||x|2 \ds' ds) ' ^ dK ^ ' 

Up to here, we didn't assume /x = constant. 
3. Necessary condition 

Now, we consider a motion with /i = constant. By the Saari's relation (jlSp . we have 

with constant v > 0. 

For the case v = 0, dx/ds — then d'^x/ds^ — 0. The equation of motion (fT6| 
yields dfi/dx = 0. Namely, the motion is homographic and the system stays one of 
the central configurations. 

Let us examine the case w > 0. In this case, the point x(s) moves on the curve 
/L((x) with finite speed v. Since the number of points dfx/dx ~ are five, we can 
always take a finite arc on which d^/dx ^ 0. To keep satisfy dji/ds = 0, the velocity 
dx/ds must be orthogonal to 9/i/9x, so we have 

(be ^ _€V_ / _djj. dfi\ 

ds \dfi/dx\ \ dy' dx) ' ^ ' 

Here, e = ±1 determines the direction of the motion. Then, the acceleration (|16p is 
given by 

_ ( £^ ( 8eu d^i\ rj \ dfi 

d? " l^(i + 4|x|2/3)|aAi/9x| \ 3\dfi/dx\'^"d^) Jd^'^ ' 



Thus, the velocity (|18p and the acceleration (|T9|) determine the curvature of this orbit 

dfi 



1 / 4C 8e / dfi\\ 3e%/7 



1 + 4|x|2/3 V V 3\dfi/dx\ \ dx 
On the other hand, the curve /i(x) — constant has its own curvature. 



dx 



(20) 




2^^^+ [If] l-f]. (21) 




\dfi/dx\'^ \\dy J dx^ dx dy dxdy 
Equate the two expressions for k, we have a necessary condition for the motion. 



AeCv 8v^ f dn 



3(l + 4|x|2/3)|9^/ax| 9(l + 4|x|2/3)|c)^/ax|2 V d 



JX 

ffdfiVd^fi ^dfidfi d^fi ^ /c)A^y6>2^\ 
3\dfi/dx\^ \\dy J dx"^ dx dy dxdy \dx J dy'^ j 

This is the condition that any motion with /i = constant and dx/dt =/= must satisfy. 
The equation of motion is invariant under the scale transformation qk — > Xqk and 
t — >■ X^/^t. This transformation makes ^/l — > Xy/l, C — > X^^'^C, x — > x, s — > X^^^'^s, 
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and V — X^/'^v. Using this invariance, we can take v = y/S without loosing generahty. 
We write C for eC. Then, the necessary condition is 



V7: 



4C 8 / 9m 



X 



V3(l + 4|x|2/3)|a/Vax| 3(l + 4|x|2/3)|9/i/5x|2V Sx 



1 / / d'^fi 



\dfj,/dyi\^ \\dy J dx^ dx dy dxdy 
and the energy is given by 




(23) 



2 



1 / d^/I\ C 



2y dt j 21 ^ ^ ' 

Substituting dy/l/dt = {dVl/dx.) ■ {dK/ds){ds/dt), dx/ds in (HH) and the 
condition (1^^ into this expression for the energy, we wih obtain the necessary 
condition expressed only by the shape variable x. However, the condition (|23p in x 
turns out to be so complex to treat. In the next section, we will rewrite the condition 
(1231) in a concise form. 



4. Invariance of the necessary condition 

Since we are considering equal mass case, the theory is invariant under the 
permutations of positions {qi}- The exchange of qi and q2 makes C — C ^^'^ ^ ~^ 
The invariance of the necessary condition ([23l) is manifest. On the other hand, the 
cyclic permutation qi q2 ^ qs ^ qi makes 

C^C'.|^.12^^ ,25) 

2 53 - 52 2 1/2 - C 



The invariance of (j23p under this transformation is not manifest. In this section, we 
will rewrite the necessary condition in a manifestly invariant form. 

4.1. Invariants 

Under the map the Lagrange points ( — ±i^/3/2 are fixed and the Euler 

points C = —3/2,0,3/2 are cyclically permuted. Let us define /ifc — I^^^/rij for 
k) = (1, 2, 3), (2, 3, 1) and (3, 1, 2). Expressions by C are. 



These three /ifc are also cyclically permuted by (|25p . Note that the exchange qi O qj 
makes the exchange /i^ /ij. Therefore, fJ. = i-ii + ^2 + A*3 is invariant under the 
permutations of qi . 

The kinetic energy for the shape change must be invariant. Actually, we can 
easily check the invariance of 

3(1 + ||C|2)2- ^ ' 

So, it is natural to treat the space of C as a metric space whose distance is given 
by the equation (P71) . and the map (PS)) is the isometric transformation. Actually 
Wu-Yi Hsiang and Eldar Straume [H [S], Alain Chenciner and R. Montgomery 



Saari's homographic conjecture 



6 



[T], R. Montgomery and R. Meckel [J showed that this metric space is 
the "shape sphere" and the distance (P7|) is the distance on the shape sphere. 
Kenji Hiro Kuwabara and Kiyotaka Tanikawa also noticed that the shape sphere is 
useful to investigate the equal- mass free-fall problem [HI [12] ■ The map makes the 
shape sphere 27r/3 rotation around the axis that connects the two Lagrange points. 
The map Q — > — C makes tt rotation around the axis that connects one of the Euler 
point (corresponds to x = 0) and one of two-body collision (corresponds to x = oo). 

Let us use the notations in the tensor analysis. We write C, — ^ ix'^, 
X = {x,y) = (x^,a;^) and di — d/dx'^. The metric tensor gij and its inverse are 



9ij 



where (5, 



3(l + ||xP)2 
J*-* are the Kronecker's delta, 

1 for i = j 







for i j. 



Let l^l be the determinant of gij, 
\g\ = det{g,j) 



16 



9 (1-f ||xP) 



As mentioned above, the configurational measure fj, is invariant, 
invariant is the magnitude of the gradient vector of fi. We write 



|VMl' = 5]5^^(a,M)(a,M)--(i + -|x 



ax 



S'^, (28) 

(29) 
(30) 

One obvious 
(31) 



Therefore, the first term of the right hand side of the necessary condition (j23p is simply 
—2C/\\7fi\. The other obvious invariant is the Laplacian of /i. 



E 



( g'- 



g\djfi 



Now, let us consider the following invariant. 



5'^"(a,/i)(a,|v/x|^ 



5/i 
9x 



d 9/i 
9x dx 



(32) 



(33) 



Explicitly performing the differentials, it yields 

3 



A = 3 1 



dx 



dx 



9_ 
16 



XI 



dfi d 
dx dx 



dfi 



d> 



Using this expression, the second and the third terms in the necessary condition ([23 



is simply expressed as, A/(2|V^|'') — A/i/|V^p. 
expressed in the following invariant form, 

r-_ 2C_ A 
The last obvious invariant what we will use is 



Thus, the necessary condition is 



(34) 



3' ' 



— A — 35) 
dx dx 
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for any invariant (j). Where, e*^ is the Levi-Civita's anti-symmetric symbol, 
r 1 for(i,j) = (l,2), 

e'' = \ I for(z,j) = (2,l), (36) 
[ for i = j- 

Then, using equations (jH)) . (IT^ and we have 
d0 e 



(37) 



4-2. Invariant variables 



For the Newton potential, it is natural to use the variables ri and r2 defined by ([S]). 
Relations between fik defined in (j26|) and ri, r2 are 

Ml = f'rVs, M2 = r^^fis, M3 = y'(f + +rl)/3. (38) 

Now, consider the expression for the above invariants |V/ip, A/i, A in terms of ri 
and r2. Let us write one of them 'ip{ri,r2). It is composed by differentials of fi by ri 
or r2 and products of ri and r2. Then, the result is composed of terms of rational 
function of y^T+rf+rf)/3, ri and r2, namely /is, /isZ/ii and /i3//i2- Then, has 
the following form 



V' = /('^i, '^2) + g(ri, r2) V 1 ^ = / — , — + 5 — , — M3. (39) 

3 M2/ VMi M2/ 

Here, / and g represent some rational functions. The function ip is invariant under 
the permutation of qi, namely the permutation of fii. So, it must be a ratio of some 
symmetric polynomials of fii. Therefore, it must have the following expression 

ip = h{p.,i^,p), (40) 

where /i is a rational function of elementary symmetric polynomials 

/i = Ml+M2+M3, = AtlM2 + M2M3 + MaMl: /5==MlM2M3- (41) 

Expression in terms of fik or in terms of /i, z^, p is not unique, since, by the relation 
psp. there is an identity fi^'^ + ^i^"^ + — 3. Namely, 

M?M2 + M2M3 + mIm? = 3/i?M2M3- (42) 
Therefore, we can eliminate using 



V = V2m7+V- (43) 



The expression oi ij; — h[p, \/l[ip + 3/9^, p) is unique. Thus, the necessary condition 
will be expressed by a function of invariant shape variables and p. 
Let us express |V/ip by p and p. In terms of r^, it is 



|VH^= ^^^"^^ +^ ^ -1^-2^^ ^^) (44) 
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^rfr^iri + ra) - rir2{rl + rl) - rfrj - ArfrKri + ra) 
+ {2rl + r\r2 - r\rl - A,r\rl - r\r\ + nr^ + 2r%) 

-^rrr^irl + rl) + 2{r\ + r^) - nra^ . (45) 

Substituting r\ — [izj [i\ and r2 = ^3/^27 we obtain, 

IV7 |2 (A'lA'i + + M3M1) / / 7 7 I 7 7 I 7 7\ 

IVmI = 6 6 6 - (MiA*2 +M2A*3 + M3Mi) 

yA^iM2A^3 \ 

-AAA^iA + A*2 + M3) - VtM2A*3(A*lM2 + M2Ai3 + M3Ail) 

+2(/i«Ai^M3 + ■••) + (AilAi2/4 + • • 0^ (46) 

In the last line, dots in parentheses represent similar 5 terms of permutation of 
Mi,/^2,/^3- Then expressing by /i, v, p, we obtain a expression. 



9p6 



As mentioned above, the expressions (|46|) and (|47l) are not unique due to the identity 
(|42p. Eliminating z/, we finally get the following unique expression 



22/i2jy3p2 _ 3j,4^2 _ 4^5p3 ^ 24/i3j/p^ + ISfiv^p^ - 27p^p'^ . (47) 



iV^il^ = -M^ + 2// + 6/ip - V - 3(2^2 _ + 3p^) V2Aip + V- (48) 

Thus, we get the expression for |V/ip in manifestly invariant variables p. and p. 
By a similar way, A/i in (ri,r2) and (/i, p) are 

(l + r2+r2)2/l d f dp\ I d f dp\ r^+rl-1 d^p 



ri dri \ dri I r2 dr2 \ dr2 I rir2 dridr- 



Ap^ p + 2p^ + 6p- 6p^/2pp + 3p'^. (49) 
Similarly, the expressions for A are 

^ ^ {l+rl+rff / dp__d_ ^ dp__d_ ^ rj + rj - 1 / dp__d_ ^ dp__d_\\ ^ 

3 \dri dri dr2 dr2 2rir2 \dri dr2 dr2 dri J J 

X = l ^4/^3 _ 24^5 ^ 32^7 _ 72^2^ + 660//p + 324/1^^ 

+36mV^ - 432p3 + + 2349Aip'' - 243p^ 

+3(24^3 _ QQ^5 _ i5g^2^ ^ 28/iV + 324^p2 



-93/iV^ - 216p3 - 27p'^p^ + \/%w + 3f^ I . (50) 
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Finally, {Dp)^ is also invariant under the exchange of g^, therefore, it has an expression 
by ^J, and p, 

^Dpf = m:!!:?^ (2(r? + rl) {rj - rlf -l) dp ^^^ dp 



SGr^rl ^''a dr2 dri 

I^Dpf = ^ ' I - (2/i + 3p)(4/ + 134/ip - 12/i3p _ 177^2 g^-^p"^) 



,2 _ pH2p + 3p) 



'2(28// + 108p - 36/iV - 45/ip2 + Mp^)y/2pp + 3p^ . (51) 



5. Proof of the Saari's conjecture 

In the previous section, we find the expression for the necessary condition (|34l) in 
terms of p and p by (|48)) . (l49t and (|50p . Since, we are assuming p — constant, time 
dependent variable is only p. Therefore, d^/l/dt — [d\/l / dp){dp/ dt). Using ([57]) . 

(d^l\\ 1 {Dpf (dvjy 

[dt J P IV/.P \dp ) ■ 

Substituting this expression and the necessary condition (j34]) into the expression of 
the energy ([M)) . we obtain the necessary condition for p with three parameters i?, C 
and /X, 

If there is some finite motion with p = constant and non-homographic, this 
condition must be satisfied by some finite range of p. However, since the right hand 
side of ([52]) is analytic function of p, the condition ([52| must be satisfied for all range 
of p. 

In the vicinity of /? = 0, we have the expansion of (j34p 



V7=ao + ai/2VP + aiP + 0(p^/2), (53) 



with 



ao 



ai/2 = 



2(l-/i^ + Cv/-l + 2/x2) 

Ml- V) 

3\/2((-2 + /i2)^-l + 2/i2- 2C(-1 + 2/^2)) 



(l-2/.2)2^/i(-l+2/i2) 



_ 3((-2 + p^){l + 6/2^) - 2C(1 + 7p^)y/-l + 2/j2) 

~ /|2(-l + 2/i2)3 ^ 



and 



|V/i|2 = /i2(-l + 2p^) - 6^/2/^5/2 ^ ^ 0(/?3/2), 

(i?/,)2 = -4/.V + 0(p'/^). 
Then we obtain the power series expansion of ([52]) by ^/p up to the order p ai p = Q. 
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The term of order p° in (j52p determine E. Therefore, this order gives no 
information for C and ^. The coefficient of order is 

-3/i5/2 (i ^ C^-\ + 2nA ' f (-2 + Ai') - 2Cv/^TTvl 
= ^ , ^. (54) 

4V2( - 1 + - Cv/^TTv) 

The solutions C of this equation are, 

C^- ' , ^£±^. (55) 
v/^TTV 2^-1 + 2^2 



For the case C = -1/V-l + V, 

^ 2^ 3^/2^3/2 9(l + 2/i2) 3/2^ , 

^TTV + (-1 + 2,2)2 (-i + 2;.2)> + Q(^^^^)- (56) 

and the order coefficient in the equation (|52|) is 

16(-l + 2,2)- (57) 



While the right hand side is always negative since , = ^J{1 + r\ + ?'i)/3 (1 + 1/ri + 
l/'''2) > 3. For the case C = (—2 + p?)/{2^J—\ + 2/i2), the coefficient 0^/2 vanish, 

and the coefficient of order in the equation ([52]) is 
3M-2 + M^) 

o-ipiTv)- ^ ^ 

While the right hand side is always positive for > 3. 

Thus, there is no parameters C and /x that satisfies the necessary condition 
This completes the proof for the Saari's homographic conjecture. 



6. Discussions 



We have proved the Saari's conjecture for equal-mass planar three-body problem under 
the Newton gravity. 

The symmetry under the permutation of the positions {qi , 52, Q3} has a clucial role 
for our method. For equal mass and Newton potential case, the necessary condition 
(|34p is a symmetric rational function of yiti, /i2 and p^. Thus, it is a function of p and 
p as in equation ([iU)) . This makes our proof simple. 

The next step will be the case with general mass ratio and general homogeneous 
potential U = X] ct > Q. For this case, however, a invariant function 

under the permutation for suffix of bodies will not have a simple form of manifestly 
invariant variables such as p and p. We hope, someday, someone may find a proof for 
the conjecture for general mass ratio under the Newton potential in some extension of 
our method. On the other hand, we are afraid that it is hard to extend our method 
to general a. We would have to find a completely new method for general a. 
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